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Abstract 



We use a one-dimensional (Id) core-softened potential to develop a phys- 
ical picture for some of the anomalies present in liquid water. The core- 
softened potential mimics the effect of hydrogen bonding. The interest in the 
Id system stems from the facts that closed-form results are possible and that 
the qualitative behavior in Id is reproduced in the liquid phase for higher 
dimensions. We discuss the relation between the shape of the potential and 
the density anomaly, and we study the entropy anomaly resulting from the 
density anomaly. We find that certain forms of the two-step square well po- 
tential lead to the existence at T = of a low-density phase favored at low 
pressures and of a high-density phase favored at high pressures, and to the 
appearance of a point C at a positive pressure, which is the analog of the 
T = "critical point" in the Id Ising model. The existence of point C' leads 
to anomalous behavior of the isothermal compressibility Kt and the isobaric 
specific heat Cp. 
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I. INTRODUCTION 



Water, the most common liquid on earth, is uncommon in many of its properties. For 
example, water at ambient pressure has an anomalous density behavior, as for T < 4°C its 
density decreases upon cooling. The subject of liquid anomalies is not limited to the density 
anomaly. Two other anomalies are the increase of isothermal compressibility Kt (density 
fluctuations) and of the isobaric specific heat Cp (entropy fluctuations) upon cooling. All 
these anomalies occur in liquid water [|^ and some occur also in other liquids [|l|-0]. 

In their pioneering work, Stell and Hemmer investigated potentials that have a region 
of negative curvature in their repulsive core ^ (henceforth referred to as core-softened 
potentials) in relation to the possibility of a new critical point in addition to the normal 
liquid-gas critical point. They also pointed out that for the one dimensional {Id) model with 
a long range attractive tail, the isobaric thermal expansion coefficient, ap = {dV/dT)p 
(where V, T and P are the volume, temperature and pressure) can take an anomalous neg- 
ative value. Debenedetti et al., using general thermodynamic arguments, confirmed that a 
"softened core" can lead to ap < [§. Stillinger and collaborators found ap < for a 3d 
system of particles interacting by purely repulsive interactions — a Gaussian potential 0. 

In this work, we investigate the relation of a core-softened potential to the liquid anoma- 
lies mentioned above. Ab initio calculations [0 and inversion of structure factor data 
revealed that a core-softened potential can be considered a realistic first order approxima- 
tion for the interaction of many materials with anomalous behavior |^,^, even in the case 
of network forming anomalous liquids [|l|. A recent work has showed that the anomalous 
behavior of a Id model can be reproduced in two dimensions {2d) as well [|T0| , |Tl]] . 

Here we will provide the details necessary to understand the preliminary results an- 
nounced in UTOf. Specifically, we investigate thoroughly an exactly solvable Id model in 



order to develop an intuitive picture of how the core-softened potential can lead to all three 
anomalies and to relate the occurrence of these anomalies to the shape and parameters of the 
potential. We also discuss the possible existence of a second critical point and its relation 
to the parameters of the model. 

The core-softened potential that we study is (Fig. |T]) 

oo < r < a 
, ^ I — Ae a < r < b 

^ ^ ' -e 6 < r < c ^ ^ 

r > c 

where r is the particle separation. The potential is composed of a "hard core" of diameter 
a which has a repulsive shoulder (henceforth referred to as "softened-core" ) of width {b — a) 
and depth Ae, and an attractive well of width (c — 6) and depth e. Unless specified otherwise 
our numerical calculations are for the choice of values a = 1, b = 1.4, c = 1.7, e = 2, and 
A = 0.5 |jl2|. A similar potential was introduced by Ben-Naim to model water anomalies 



A continuous potential of an analogous form models the interaction between clusters 
of strongly bonded pentamers of water |[T^. This type of interaction qualitatively repro- 
duces the density anomaly of water. At sufficiently low pressures and temperatures, nearest 
neighbor pairs are separated by a distance r 6 to minimize the energy. As temperature 



2 



increases the system explores a larger portion of the configurational space in order to gain 
more entropy. This includes the penetration of particles into the softened core, which can 
cause an anomalous contraction upon heating. 

The relation between the density anomaly and the shape of the potential can also be 
understood using the thermodynamic relation 

{dV/dT)p = (3^ {P{{5Vf)pT + {5V5E)pt) (2) 

where (3 = l/ksT with kp being the Boltzmann factor, E is the energy and {■■■)pt is the 
thermodynamic average in a constant P, constant T ensemble. For a system with density 
anomaly the left hand side of Eq. (Q) is negative at temperatures below the temperature of 
maximum density Tmd- The first term on the right-hand side, proportional to the square 
of the volume fluctuations, is always positive. Thus below Tmd the second term must be 
negative, which requires anti-correlation between the fluctuations in E and V. This kind 
of anti-correlation exists for the core-softened model when the fluctuations occur between 
the states with r ~ 6 and the states with r ^ a which have a larger energy but a smaller 
volume. We conclude that the core-softened potential is a candidate for generating the 
density anomaly. 

The paper is organized as follows: In Sec. |I| we present the exact solution for the 
Gibbs potential. Sections III-VI discuss the anomalies in the density and entropy, and 
their response functions compressibility and specific heat. In Sec. VII we discuss interesting 
analogies with the Ising model. Finally, Sec. VIII interprets the anomalies in terms of two 
different local structures. 



II. EXACT SOLUTION FOR GIBBS POTENTIAL 

To derive the Gibbs potential for the model, we choose c < 2a to restrict the interactions 
to nearest neighbors. The Id model is then exactly solvable using standard methods 0,|T^, 
16[. The partition function is (Appendix ^ 



where N is the number of particles, AV is a discretization factor which is needed for keep- 
ing Z dimensionless and does not enter the equation of state and A(T) is the De Broglie 
wavelength. Here 

Q{T, P)= dx e-^^"e-^"(") (4) 
Jo 

where u{x) is the interaction potential. 

The Gibbs potential G{T,P) = —ksT In Z{T, P) is then 

G(T, P) = -kpT [(A^ - 1) In ^](^, P) - In A(T) - 2 In P - In AV] . (5) 

In the thermodynamic limit oo, the terms 21nP and InAV^ are negligible compared 

to the other extensive (order A^) terms. 
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Equation is valid for any Id system in which each particle interacts only with its 
nearest neighbors. If the interaction potential is given by Eq. (|l|), then 

fi(T,P) = ^vI;(T,P) (6) 

where (Appendix ^ 

*(T, P) = [<f)^ea + (0 - + (1 - m) (7) 

and e^{P,T) = e^^^^ (f){T) = 

III. DENSITY ANOMALY 

We calculate the equation of state using the definition V = dG(T, P)/dP. From Eq. @ 
we find that, in the thermodynamic limit N oo, the average "volume" (length in Id) per 
particle is 

' n{T,p) dP 

Using Eqs. @ and (H) we find the equation of state 

_kBT _ ksT d^{T,P) 

^ ~ ~T~ ~ m{T,p) ap 

where 

^^^^ = -(3 (cjy^ae^ + (0 - <p^)he, + (1 - <p)ce,) . (lo) 

In the high temperature limit where — > 1 and 6* — > 1, the equation of state (^ tends to 



(9) 



P 



+ a, (T^oo) (11) 



the equation of state for an "ideal gas" of rods, that is a system of non-interacting rods of 
length a. 

At T = 0, w as a function of P has a discontinuity at 

P-P =(lll^ (121 

which we call the upper transition pressure. At T = 0, 

v = l' Tp"' (13) 

la P> Pup ^ ^ 



We will return to this first-order transition in Sec. |VI1I| . 
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Next we study the equation of state for fixed pressure. For each value of P, we find the 
value of v{T) using Eq. @ (Fig. |^). The isobars separate into two different groups: 

(i) For P > Pup; V increases monotonically with T, starting from its minimal value 
V = a. At higher pressures, the nearest neighbors are pushed inside the softened core and, 

result, the density anomaly does not occur. 

(ii) When P is lowered just below Pup, f = 6 at T = [Eq. (0)], and the v{T) isobars 
show a maximum at a temperature of minimum density T^d and a minimum at Tmd I^j^], 
T3|, |14|, |17| . For P < P^p, the system starts at T = from the bottom of the energy well at 



V = b. Upon heating, the particles first start to explore the wider region in the bottom of the 
potential well and the system expands. At higher temperatures, the particles penetrate the 
softened core and the system shrinks, showing an anomalous temperature-driven contraction. 
For even higher temperatures, the particles move outside of the potential well, causing the 
system to expand. 

The minimum density point T^d is of interest. In real liquids with Tmd, T^d is rarely 
observed possibly because it would be located at a very low temperature, where the 
liquid phase is not stable. As pressure is lowered further, the maximum and minimum 
density points coincide at some point (Tiow,Piow), which can be found from the system 
of equations {dv/dT)p = and {d'^v/dT^)p = 0. We have observed this behavior upon 
changing the parameters of the interaction potential (see. e.g.. Fig. ^). For P < Plow, no 
density anomaly is observed, as shown in Fig. ||b. 



IV. ISOTHERMAL COMPRESSIBILITY ANOMALY 

The isothermal compressibility is defined as 

Kt is thus the response of the volume to its conjugate variable pressure and it is proportional 
to fluctuations in specific volume 

Kt oc {{SVf). (15) 

For most materials {dKx / dT)p > 0, so fluctuations decrease upon cooling. In the 
case of water, for a wide range of pressures, Kt passes through a minimum and shows an 
anomalous increase upon cooling. Along the P = latm isobar, for example, the minimum 
compressibility point is around 46° C |Q. 

From Eq. (|^) and Eq. (p!4D, we find 

K - ^ ( ^ I ^ d'nT,P) ( 1 dnT,P) y\ 

^ Pv\p^^{T,p) ap2 l^(T,p) dP I r ^ ^ 



where 



^^f^ = {<p'a'ea + (0 - <P')bX + (1 - <P)c%) . (17) 
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Figure shows Kt as a function of T along isobars from Eq. (plGl). As T — > oo, Kt tends to 
its ideal gas value 1/-P, which is also predicted by Eq. ([TTl) . As T — > 0, 
using Eq. (p!^, we find 

Kt^^ (T-^0,P^P,p). (18) 

As T decreases from infinity to zero, Kt passes through two extrema, a minimum and a 
maximum, between which {dKx / dT)p < 0. Thus the Id core-softened model generates a 
compressibility anomaly. 



V. DENSITY AND COMPRESSIBILITY EXTREMA LINES 



The locus of the points Tmd in the P — T phase diagram is of special interest. For water, 
the shape of the Tmd line helps distinguish between different scenarios proposed to explain 
water's anomalies |18|^0|. In simulations ||2l|-^, the Tmd line presents a "nose," i.e., a 



point in which as P decreases the slope changes from negative to positive passing through 
an infinite value (Fig. 1). 

In Fig. |, we present the Tmd and T^d with the same parameters as in Fig. ||. We 
observe that both lines originate at Pup p5| and terminate at Piow Moreover, the Tmd line 
has a negative slope for large P which is in agreement with experimental water and for lower 
P the slope of the Tmd line changes sign and so a nose is present. 

As shown in Fig. H, the locus of extrema in compressibility intersects the Tmd line at its 
nose, as predicted by Sastry et al. using thermodynamic relations. Moreover wherever 
the Tmd line has a negative slope, the compressibility must increase upon cooling in the 
region to the left of the Tmd line |T^, resulting in a line of Kt maxima which originates 
from the point C'iT = 0, P = P^p). 

Noteworthy is the existence of a starting and ending point for the Tmd line. The starting 
point is the point C'{T = 0, P = Pup), and the ending point (Tiow, Plow) where T^d and Tmd 
meet. In Fig. we show that the overall phase diagram does not change qualitatively upon 
varying the parameters of the potential. 



VI. ENTROPY ANOMALY 



Since {dS/dP)T = -{dV/dT)p, if {dV/dT)p < 0, then 

> 0. (19) 



dS_ 
dP 



T 



Equation ( |T^ ) is anomalous because, contrary to intuition, compressing the system at 
constant T increases its entropy [^]. To calculate the entropy from Eq. (^, we use 
S = -{dG/dT)p 
We obtain 

^ +H^iT,P)/iApP))--J—^^%^ (20) 



Nks 2 V V ' //V A- ^(T,P) df3 
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where 



= (eA - Pa)ct>''ea + e(0 - \ct>^)e, - (0 - 0^)P6^, - ecj^O, - (1 - ct>)Pce,. (21) 

Using Eq. (|20|) , we plot the entropy for two isotherms (Fig. ^). For T = 0.6, there is no 
density anomaly (Fig. |^a) and hence no entropy anomaly. For T = 0.5, there is a density 
anomaly and hence an entropy anomaly; Fig. shows this anomalous increase of as a 
function of P. 



VII. ANALOGIES WITH THE ISING MODEL 

A. First Critical Point 

The gas-liquid first order transition line ending at a critical point which is present in 
higher dimensions shrinks to the point C = (T = 0, P = 0) for a Id fluid in which the 
particles interact with an attractive potential. This point C is the remnant of what is a 
critical point in higher dimensions: for example Kt diverges for T — > 0, 

i^T-l/T (r^O,P = 0), (22) 

analogous to the divergence of the magnetic susceptibility xt along the zero field line H = 
for the Id Ising model, 

Xt~1/T {T^0,H = 0). (23) 



The constant-pressure specific heat Cp = T{dS/dT)p is obtained from Eq. (P0|), with the 
result 

C^^S J 1 d'^{T,P) _ ( 1 d^{T,P) y\ 
Nkn 2 ^ \^(T,P) V^(T,P) d(3 J j ' ^ ^ 

For small P and T, the specific heat has the form 

Cp ~ ao + a2{f3Pf (25) 

the analog of the Ising case for small H and T p7[| 



Ch ao + a2i(3Hf. (26) 
These features are common to model Id fluids with attractive potentials. 

B. Second Critical Point 

Next we will discuss the "remnant" of a second critical point for our core-softened model. 
In addition to the divergence along the P = isobar, we find a divergence along the P = Pup 
isobar (Fig. 3), 
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[T^O,P = Pg. (27) 

We also observe from Eq. (^) and Fig. |^ that there is a discontinuity in the order 
parameter v when crossing C'{P = Pup, 7" = 0) along the T = axis; this is the analog of 
what happens to the magnetization when crossing the H = point along the T = axis for 
the Ising model. 

Next we consider Cp. Taking the limit T ^ of Eq. (^) and defining H = \P — Pup|, 
we find 

where we have introduced the parameters A = b — a and w = exp{—i3AH). Eq. ( PBD is 
the analog of Ch for the Ising model. Figure ^ shows the anomalous behavior of Cp as a 
function of T along different isobars. 

It is interesting to note that simulations of water using the ST2 potential display a 
compressibility anomaly due to the presence of a second critical point in the metastable 
region of the liquid p2|,p8|, and a line of Kt maxima originates from this second critical 
point. 



VIII. INTERPRETATION IN TERMS OF TWO LOCAL STRUCTURES 



The density and compressibility anomalies can be related to the interplay between two 
local structures: an open structure in which the nearest neighbor particles are typically at a 
distance b, and a denser structure in which the nearest neighbors penetrate into the softened 
core and are typically at a smaller distance a. The favored local structure is determined by 
the Gibbs potential per particle 

g{T, P) = mm{u + Pv-Ts}, (29) 

which is shown in Fig. ^ as a function of the volume per particle, at T = for two different 
values of P. For P < Pup and T = 0, the minimum corresponds the "open structure" with 
r ^ b. 

Increasing P increases the value of -u + Pv — Ts for the open structure (Fig. |^). For 
P > Pup, the minimum corresponds to the "dense structure" with r ^ a. 

The value of Pup of Eq. ( [T^ ) can also be found by equating u + Pv — Ts (at T = 0) for 
the two local minima, which results in 

Pup (^opcn ^donso) / (^opon ^donso)' (^^) 

Substituting Wopen = — e, ^dense = " Ae, fopen = b, Vdense = CL rcsults in the same expression as 
Eq. (pD. 

For higher dimensions, Eq. (pOf ) helps to estimate the pressure region in which a transition 
between a dense and an open structure could happen. 

For d = 1, the contribution of the Ts term makes the double well structure of Fig. 4 
disappear when T > 0. This may not be true in higher dimensions. If we assume that the 
qualitative shape for d> 1 changes little from the d = 1 case, then for > 1 there can exist 
a first-order transition line for small T, eventually terminating in a critical point C . 
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IX. SUMMARY 



We used a Id core-softened potential, which mimics the effect of hydrogen bonding, to 
develop a physical picture for some of the anomalies present in liquid water. We discussed 
the relation between the shape of the potential and the anomalies in density and entropy 
and their associated response functions Kt and Cp. The form of the potential leads to the 
existence at T = of a low-density phase (favored at low pressures) and of a high-density 
phase (favored at high pressures), and to the appearance of a point C at a positive pressure, 
which is the analog of the T = "critical point" in the Id Ising model. 
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APPENDIX A: DERIVATION OF THE FREE ENERGY FOR ID 
CORE-SOFTENED MODEL 

We start from the general definition of the partition function 

1 c 
Z{T, P) = ^^^^^,^jv J dVe-f^^^ J d^pd^x e-/5W(piv..,Piv,xi,...,x^) (^1) 

where N is the number of particles, V is the Id system size, Xi and pi are the position 
and momentum of particle i respectively, Hi^pi, ...,pn,xi, ...jXn) is the Hamiltonian of the 
system. AV is a discretization factor for V, which is needed for keeping Z dimensionless and 
does not enter any equation of state. The Hamiltonian can be partitioned into its kinetic 
and potential parts. 

H{pi, ...,pn,xi, ...,xn) = XI ^ + U{xi, ...,xn) (A2) 

i=i 



One can separate Eq. (|A1|) into a momentum and a configurational integrals, where the 
momentum part is 

Z,(T).^/d^exp(-aE||). (A3) 

This integral can be written as the product of N integrals over momenta and then using 



Gaussian integral formula p9| we find 

Z,{T) = ^ dp exp (-/3pV2m) Y = [A(T)]-^ (A4) 

where A(T), which has the dimension of length, is the temperature dependent De Broglie 
wavelength A(T) = h/ \/2TimkBT . Thus the partition function takes the form 
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In order to take further steps, we use the fact that the range of interaction is less that 
twice the hard core (c < 2a). As a result, one can think of the particles in the Id system 
as a chain, in which each particle interacts only with its nearest neighbors. Thus for each 
arrangement of the particles, we can write the interaction potential U as the sum of — 1 
terms 

U{xi, ...,xn) = U{xn - xm-i) + U{xn-i - XM-2) + ••• + U{x2 - xi) (A6) 
Using Eq. (^) we rewrite Eq. (^) as 

Z{T,P) = ^^y^^^rp^N Jo Jo dxN-iE{xN - xn-i) 

/ dx]\!^2E(xj^^i — xn^2) ■ ■ ■ dxiE(x2 — Xi), (A7) 
Jo Jo 

where E is defined as 

E{x) = e-^^^'^K (A8) 



In Eq. ([A7|) we have used the symmetry of Z under permutation of the particles and thus 
have written Eq. (|A7|) only for a specific order of particles, in which xi < X2 < ■ ■ ■ < xn- 
All other permutations are equivalent to this integral, after renumbering the particles. This 
results in a permutation factor A^! which has canceled the A^! factor in the denominator. We 
rewrite Eq. ( |A7| ) as 



{AV)A{Ty 

where we have introduced the function Ei{x), defined recursively as 

F^{V) = r dxj,FN^,{xN) = (1 * Fn-iW) (AlO) 
Jo 

Fi{xi+i) = / dxiE{xi+i - Xi)Fi_i{xi) = {E* Fi^i){xi+i) (All) 
Jo 

Foixi) = 1. (A12) 
The star operator *, represents the binary convolution functional defined as 

{f*g){x)^ r dyf{x-y)g{y). (A13) 
Jo 

We further use the notation £[/] for the Laplace transformation functional which is defined 

as 
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oo 



£[f]{z) ^ e-^^f{x) (A14) 







to find the following simple form for the partition function 



1 

{AVjA{Ty 



z{T,P) = twttwt^^IFnKPP). (A15) 



Note that according to Eqs. ([XToD , (^T|) and ([ATsD 



N-l 



= JKvjMTj^^^^ *{E*{E*...{E*1).. •))](/?P). (A16) 
Next we use the convolution theorem, which states that the Laplace transform of the con- 



volution is equal to the product of the Laplace transforms of each function |^ 

£[f*g]{z)=C[f]iz)xC[g]iz). (A17) 
and also the formula for the Laplace transform of a constant function p3 



to obtain 



C[c]{z) = ^ (A18) 



ZiT, P) = 7WT7V777FVW77^MT. P)f-'- (A19) 



(AV)A(T)^(/3P)2 
Here fi(T, P) is defined in Eq. (§. 



APPENDIX B: LOW-T LIMIT 

In order to derive properties of the free energy at low temperature approaching the point 
C'{T = 0,P = Pup), we note that in Eq. (|^ as T — 0, for small P, the 06'f, term dominates, 
while for large P the (f)^6a dominates. In order to find the limit of \E'(T, P) and its derivatives, 
we note that 

^=exp[/3(P(6-a)-(l-A)e)]. (Bl) 

The value of P where the (pOb term balances the 4>^6a term follows an equating the argument 
of the exponential to zero in Eq. ( |B1D . The result for Pup is given in Eq. (T^). Using (|T2p 
and ( |B1| ), Eq. (|^ becomes 

^(T, P) ~ re^{l + exp {-(3{b -a)\P- P^pD) (B2) 

where 

a = 1,6^ = Ob (P < Pup) .r^o^ 

a = X,e, = ea (P>Pup) ^ ^ 
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For a^(T, P)/dP we use Eq. (|B|) and dd^/dP = -^xO^ to find 

d^{T,P) 



dP 



(/?x)vI/(r,P), (B4) 



where x = h ioi h < Pup and x = a ioi h > P^p. Using Eqs. ( |B3| ) , (p^) and (p^ ) in Eq. 
leads to Eq. (|l|). 

To find the scaling form of Kt for T ^ 0, we use the equation for the second derivative 
of ^(T, P) 



d^^{T, P) 



(x{l3xy^{T,P). (B5) 



Qp2 

^From Eq. ([T6|), we find Kt for any P 7^ P^p to behave as 

i^T ~ iPvP'')-' = {P + Px P^r' (T ^ 0, P ^ Pup) (B6) 

which is consistent with Eq. (p!8|). 

In order to derive the limiting expressions for the entropy and specific heat, we must 
differentiate the free energy with respect to f3. We start by rewriting Eq. (^) as 

g{T,P) = — = keT ln{A/n) = keT (3 -\n^!{T,P) + \nP + constants] . (B7) 

For ^'(r, P) we rewrite( ^ as 

^(T,P) ~ P(T,P) X (1 + w(T,P)) T^O (B8) 

with the definitions 



P(T,P) = exp {^{ta-Px)) 

wiT, P) = exp(-/3)(6 - a)|P - P^pl 



(B9) 



Using these definitions we find 



a^(T, P)/dp oc {ea - Px)E (1 + e(l - T)) 
a2^(T, P)/a/32 oc (ea - Pa;)^^ (1 + e(l - T)^) 



(BIO) 



^ ^ - Pxf (i±^ - ^i^^^) (BID 

2 V 1 + w 1 + wY I 



Using the above equations and Eq. (|2^), we find 



whose leading term is Eq. 
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FIGURES 



FIG. 1. General form for the core-softened potential studied here. The length parameters 
a, b, c and energy parameters e, A are shown. 

FIG. 2. Isobars of v, the average length per particle, for the discrete Id core-softened po- 
tential showing the conditions under which a Tmd and TmD exist, (a) The parameter values are 
a = l,b = 1.4, c = 1.7, e = 2 and A = 0.5. From Eq. (|l2|) , these values result in Pup = 2.5. Piow is 
almost zero, so no P < Piow isobar is shown. The Tmd point is marked by a filled ellipse and the 
Tmd point by an open ellipse, (b) The parameter values are a = l,b = 1.2, c = 1.7, e = 2, A = 0.5 
and thus from Eq. (|l2|), Pup = 5. Now Piow ~ 3.3, so a P < Piow isobar is shown. 

FIG. 3. Isothermal compressibility for the same parameter values as Fig. 2a. (a) Isothermal 
compressibility along different isobars, with their maxima marked by filled circles. Kt along the 
isobar Pup diverges as T — > 0. (b) Log-log plot of the same results showing the divergence of Kt 
as 1/T along the critical isobar P = Pup. 

FIG. 4. (a)The loci of the two density extrema {Tmd and TmD) and the two Kt extrema 
[K^^^ and for the discrete potential of Fig. 1 with the same parameter values as Fig. 2a, 

so Pup = 2.5 and Piow ~ 0. The locus of Kt extrema has two "loops," one connected to the point 
C'{T = 0, P = Pup) and the second to the point C{T = 0, P = 0). (b) Same for parameter values 
of Fig. 2b, for which Pup = 5 and Piow ~ 3.3. The two "loops" join. 

FIG. 5. (a) The region of the P — T plane where density and entropy are anomalous (grey); 
the parameters are the same as in Fig. 5a. (b) The behavior of the entropy along two different 
isotherms. The T = 0.5 isotherm intersects the anomalous region and shows the maximum and 
minimum marked by the circles. The T = 0.6 isotherm is outside the anomalous region and does 
not show any anomaly. 

FIG. 6. Constant-pressure specific heat along different isobars, P = 2.7 = Pup -|- 0.2 and 
P = 2.3 = Pup — 0.2. Near T = 0, the two curves coincide, as predicted by Eq. (pS]), since the 
values oi H = \P — Pup| are identical. 

FIG. 7. The function u + Pv — Ts at T = for the core-softened potential of Fig. |^. The 
equilibrium value of v{P) is determined as the absolute minimum, which is located at v = b for 
P < Pud and at f = a for P > Pun. 



15 



u(r) 



-'ke 
-e 




17 




18 




19 





dS/dP<0 


1 

1 


(a) 




dV/dT>0 


1 
1 
1 


- 






1 
1 

4 




> 




1 
1 






dS/dP>0 


1 
1 






dV/dT<0 


1 
1 

1 








1 
1 








1 
1 
1 
1 

1, 





0.0 0.2 0.4 rp 0.6 0.8 




-0.4 I ^ • ^ • ^ ' 

1 2 p 3 



20 




21 




22 



